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ABSTRACT 

Using analytic arguments and a suite of very high resolution 10'^ M0 per particle) 
cosmological hydro-dynamical simulations, we argue that high rcdshift, z ~ 10, M 
10* Mq halos, form the smallest 'baryonic building block' (BBB) for galaxy formation. 
These halos are just massive enough to efficiently form stars through atomic line 
cooling and to hold onto their gas in the presence of supernovae winds and reionisation. 
These combined effects, in particular that of the supernovae feedback, create a sharp 
transition: over the mass range 3- 10 X lO^M© , the BBBs drop two orders of magnitude 
in stellar mass. Below 2 x IO^Mq galaxies will be dark with almost no stars and no 
gas. Above this scale is the smallest unit of galaxy formation: the BBB. 

We show that the BBBs have stellar distributions which are spheroidal, of low 
rotational velocity, old and metal poor: they resemble the dwarf spheroidal galaxies 
(dSphs) of the Local Group (LG). Unlike the LG dSphs, however, they contain signifi- 
cant gas fractions. We connect these high redshift BBBs to the smallest dwarf galaxies 
observed at z = using linear theory. A small fraction (~ 100) of these gas rich BBBs 
at high redshift fall in to a galaxy the size of the Milky Way. We suggest that ten 
percent of these survive to become the observed LG dwarf galaxies at the present 
epoch. This is consistent with recent numerical estimates. Those in-falling halos on 
benign orbits which keep them far away from the Milky Way or Andromeda manage 
to retain their gas and slowly form stars - these become the smallest dwarf irregular 
galaxies; those on more severe orbits lose their gas faster than they can form stars 
and become the dwarf spheroidals. The remaining 90% of the BBBs will be accreted. 
We show that this gives a metallicity and total stellar mass consistent with the Milky 
Way old stellar halo. 
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1 INTRODUCTION 

The Local Group (LG) of galaxies provide a unique test- 
bed for galaxy formation theories and cosmology. Their close 
proximity allows individual stars to be resolved giving accu- 
rate kinematic s, stellar population s a nd star formati on his- 
tories (see e.g. iKlevna et al.ll200il and lDolp hinll2002ft : their 
spatial distribution can be compare d with cos mological pre- 
dictions to give useful constraints iMoore et al . 2006); and 
their large mass to light ratios can be used, through dynam- 
ical mo delling, to place co nstraints on the nature of dark 
matter jKlevna et al.ll200j) . 

The LG dwarf galaxies are usually split into three types: 
dSph galaxies, which typically have old stellar populations, 
are spheroidal in morphology, lie close to their host galaxy^ 



and are devoid of HI gas; dirr galaxies, which have younger 
stellar populations, irregular morphology, lie further away 
from their host galaxy and contain significant HI gas; and 
the transiti on galaxies, which are in between the dSph and 
dIrr types llMatedllQQa) . 

In our current 'vanilla' cosmological paradigm (ACDM - 
cold dark matter with a cosmological constant) , all structure 
forms from the succe ssive mergers of smaller substructures 
JWhite fc Reejll97d) . While this theory has been tremen- 
dously successful on scales larger than ~ 1 Mpc, on smaller 
scales it has fared less well (see e.g. ..D'O nghia & Lake 2004i). 
A now long-standing puzzle is the 'missing satellites' prob- 
lem: there appears to be an order of magnitude fewer 
satellite galaxies in the LG than would naively be pre- 



* Email: justin@physik.unizh.ch 

^ We use the terminology 'host galaxy' throughout this paper 



to refer to either the Milky Way (MW) or Andromeda (M31) de- 
pending on which of these is closer to the satellite being discussed. 
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diet ed from the mass fun ctio n of dark matter h alos (see 
e.g. iKauffmann et al ] |l993l and iKlvpin et al1ll999h . A num- 
ber of solutions to this puzzle have been presented, the two 
main threads bein g ei ther to alter the nature o f dark matter 
jBode et alJl200ll and lAvila-Reese et alJl200ll) . or to invoke 
some form o f feedback from sup^ernovae explosions (see e.g. 
lLarsorll974l an d Efstathiou 200 (i'), or photo-evapora tion (see 
e.g. iQuinn eral.lll996l and iBarkana fc Loebl^l999^ . In the 
feedback scenario, one might naively expect only the most 
massive substructure satellites, with the deepest potential 
wells, to form stars and remain visible in the LG at the 
present epoch JStoehr et alJl2002^ . However, this presents 
a problem since the most massive substructure dark mat- 
ter halos predicted by ACDM models have central stellar 
velocity dispersions which are factors of 2-3 larger than 
those observed in the LG satellites, oven after extreme tidal 
stripping of shocking o f t hese ha los (Havashi ct al. 2003,, 
iKazantzidis et al.ll20(3j and lRead et al.ll2006aD . 

An alternative view, is that the L G satellites are fossil 
galaxies left o ver fr o m reionisation iBuUock et alJ l200l 
BuUock et all l200ll fe enson et al.' '2002', iKravtsov et^ 
2004 ISusa fc Urnemur a 2004, Ka wata et al.. 20051 
Ricotti fc GnedinI |2005[ and iGnedin fc Kravtsovl l200d) . 
In this scenario, only those rare over-dense peaks which 
collapse before redshift z ~ 10 (the epoch of reionisation) 
and achieve a potential well deep enough to form stars 
remain visible at the present epoch; the remaining satellite 
galaxies have star formation quenched by the background 
UV flux from reionisation. In this model it is not the 
most massive substructure halos at 2: = which are the 
LG satellites, but rather the survivors from halos which 
form stars at z; ~ 10. There is observational evidence for 
such a scenario, bot h indirectly from high redshift quasar 
absorption spectra JWvithe fc Loebl l2006^ ■ and from the 
star formation histories of LG dirrs, whi ch show a stron g 
suppression in star formation up to 2; ~ 1 llSkillmanll2003) . 

While this general model has been investigated by a 
number of authors in th e literature, th e re are significant 
differences in the details. IKravtsov et alJ (|20o3) argue that 
only a few of the LG satellites are genuine fossils; the ma- 
jority were significantly more massive in the past, formed 
most of their stars after reionisation (at z ~ 3), and then 
subsequently lost their mass through tidal stripping and 
shocking. They find that the central velocity dispersions 
can be sufficiently lowered in their model once satellite- 
satellite interactions are taken into account alongside strip- 
ping and shocking from the Milky Way (MW) or Andromeda 
( M31). This model com ple ments the 'tidal mode l' proposed 
bv'Mav er et al.l ||2001a') and'Ma ver et al.l l|_2001b') to explain 
the distance-morphology relation between dSphs and dIrrs. 
In these two papers, it is suggested that all of the Local 
Group satellite galaxies started out looking more like the 
dIrrs with a disc-like morphology. Those satellites on orbits 

^ The epoch of reionisation is caused by UV flux emitted by the 
first forming massive stars. Throughout this paper we suggest 
that reionisation occured at 2: ~ 10. Observationally there are 
two bounds on t his epoch. The old WMAP satellite data favour 
larger redshifts ISoergel et alJl2003^ : th e data from Qu asar ab- 
sorption spectra favour smaller redshifts <Fan et al .l2005t but see 
also .Bunker ct al. 2006); the new WMAP data favour our chosen 
redshift iSpergel et aljl200d) . 



which brought them close to the MW or M31 then formed 
induced bars which buckled leaving a spheroidal remnant: a 
dSph^ 

iRicotti fc GnedinI and iGnedin fc Kravtsovl 

(2006) focus on the pure fossils left over from reionisation. 
They present a detailed cosmological model which includes 
radiative transfer and molecular cooling from H2. These 
new key ingredients allow them to study star formation 
in mini halos with virial temperatures T < 10'* K - the 
temperature at which hydrogen starts to become collision- 
ally ionised. They find, contrary to previous studies, that 
such mini-halos can cool efficiently and form stars. They 
posit that such halos could then (if they survive) be the 
progenitors of dSph galaxies in the LG. They require more 
massive halos to then become the dirr galaxies. However, 
they do naturally recover the spheroidal morphology, 
low gas fractions, low rotational velocity and old stellar 
popu lations observed in the LG dSphs (see e.g. iMated 
Il99^. without recours e to any tidal tr ansformations. 

[ Ricotti fc GnedinI (I2OO3) and IGnedin fc Kravtsovl 
i2006l) studv in detail the effects of photo-ionising feedback 
from star formation, but they do not include the effects of 
feedback from supernovae winds. It is important to separate 
feedbac k from supernovae (whic h is i ncluded in the mod- 
els by iRicotti fc GnedinI ll2005ll and IGnedin fc Kravtsovl 
i2006ll l from feedback from supernovae winds, which are 
not. In this paper, we implement and test the effect of both 
kinds of feedback. This is a key difference in the work we 
present here. Observationally, we can see that supernovae 
winds are driven in galaxies undergoing a phase of star 
formation and that such winds are i mportant, e specially on 
the scale of dwarf galaxies (se e e.g. l Ott e t alj|20 05 | ). The- 
oretic ally, iMac Low fc Ferraral 1119991 ) and iMarcolini et alJ 
i200(j) have shown that there is a transition scale at 
about Merit ~ 10^ — 10* Mq below which dwarf galaxies 
efficiently lo s e the ir g as from supernovae winds'^; while 
iDekel fc Sil^ Jl98dl and lDekel fc Wool ll2003ll demonstrated 
that supernovae feedback could account for the global 
scaling relations of the LG dwarfs. 

In this paper, we use a suite of very high resolution cos- 
mological hydro-dynamical simulations, which include gas 
cooling, star formation, feedback from supernovae, galac- 
tic winds and reionisation, to study a new model for the 
formation of the LG dwarf galaxies. We do not include 
the effects of the detailed radiative transfer and cooling 
physics required to model star formation in mini-halos with 
mass M < 10^ Mq, since we are interested primarily in ha- 
los more massive than this. We suggest that the smallest 
dwarf galaxies of the LG share a common progenitor: rare, 
~ 3(T, ^ 10** Mq, dark matter halos at z ~ 10. These ha- 
los are just massive enough to efficiently form stars through 
atomic line cooling and to hold onto their gas in the pres- 
ence of supernovae winds and reionisation. As a result they 
are the smallest 'baryonic building block' (BBB) available 
for galaxy formation. 

Some of these gas-rich early forming galaxies fall-in late 
to the LG and survive as dwarf galaxies. Those in-falling 



^ We note, however, that this conclusion is degenerate with the 
mechanical luminosity in the wind. For the mean plausible range, 
efficient mass loss occurs at ~ 1O^M0. 
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halos on benign orbits which keep them far away from the 
MW or M31 manage to retain their gas and slowly form 
stars - these become the dirrs; those on more severe orbits 
lose their gas faster than they can form stars and become 
the dSphs. This suggests that the dIrrs should also have an 
old extended spheroidal component of stars: a stellar halo. 
There is incre asing observational evidence that this is in- 
dee d the case ('Minnit j_fc Ziilstra ''l99ef. 'Minniti et al."2003l 
and lAparicio"et aL.200Ci and see Vaduvcs cu et al. ,2005i for 
a study of dIrrs outside of the LG). It also suggests that the 
star formation histories of the LG dwarfs should all show 
an early, pre-reionisation burst of star formation. This also 
appears to be the case (see e.g. iHernandez et al1l2000l and 
iDolphinlEoO^) . 

The idea of a sm allest building block has a long his- 
tory in the literature. iHovld (119531) pointed out that gas 
cooling becomes efficie nt for ionised hydrogen at 10"* K; 
iPeebles fc Dickj il968l) presented a model with a smallest 
mass block of lO ^Mp^ fo r star formation in the context of 
globular clusters: iPeeble s (1984!) updated this argument to 
include dark matter and showed that the relevant mass was 
~ 10* Mq, as suggested her e : and , from an observational 
point of view, ISearle fc ZinnI (Il97^ showed that the galac- 
tic globular cluster abundances are not correlated with dis- 
tance - a result which has led to the hierarchical merging 
model being currently fa voured over t he then-popular model 
of monolithic collapse iEggen et alj |1962) . Here, we take 
the next logical step by investigating the morphology and 
kinematics of the stars, gas and dark matter in these high 
redshift BBBs, and making a link to the smallest galaxies 
observed in the Local Group at the present epoch. 

To test our model, we use a small box size of 1 Mpc 
and stop the simulation at 2; = 10 (to avoid simulations 
becoming non-linear on the scale of the box). With such 
a small box size we achieve an unprecedented mass resolu- 
tion of ~ 10^ Mq per particle. This allows us to accurately 
track the kinematics and morphology of galaxies of total 
mass M ~ 10* Mq and compare our results with observa- 
tions from the LG dwarfs. Ideally, one would like to use a 
larger box and evolve all the way to redshift z = 0. How- 
ever, this is not technically feasible for a mass resolution of 
~ 10^ Mq per particle, at the present time. Instead we have 
to compromise by stopping at redshift z = 10 and making 
the link to redshift z — using a mixture of linear theory 
arguments (section P)l and re sult s from other studie s in the 
literature dMoore etal.ll2006l and lMaver et al.ll2005^ . 

This paper is organised as follows: in section |2| we 
present some theoretical motivation for our model. In sec- 
tion |21 we describe the simulations. We used a control sim- 
ulation with only dark matter, and five other simulations 
which explore the effect of the star formation prescription 
and supernovae winds of varying strength. All simulations 
were run with identical initial phase space distributions. In 
section |1| we describe the results from our suite of high res- 
olution simulations. In section |S] we discuss our results in 
the context of recent observations. Finally, in section |S| we 
present our conclusions. 
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Figure 1. The number of halos, N, formed up to redshift z = 
15, 10, 5 (solid, dotted and dashed lines), in a given mass range 
which are likely to have fallen into a larger halo of mass Mq = 
2 X 10^-^ Mq at zq = 0. Over-plotted are the number of dSph 
galaxies currently observed around the MW (horizontal line) and 
the mass scale M^^^ (vertical line). 



2 THEORETICAL MOTIVATION: THE 
'BARYONIC BUILDING BLOCK' 

The mass scale, Afcrit ~ 10* Mq , is interesting. Firstly, 
it corresponds to the mass below which a galaxy can- 
not retain its gas in the presence of supernovae winds 
(Mac Low fc Fe rrara 1999; Marcolini ct al. 2006*). Sec- 
ondly, the viria l temperature of halos at a given redshift, 
Tv, is given by (|Peacocklll999l) : 

n/K = io^-\M/io^^MQf/\fcnh'f^'\i + z) (1) 

where M is the mass of the halo; z is the redshift; h = 0.72 
is the dimensionless Hubble parameter at the present epoch; 
Q — Q{z) is the ratio of the halo density to the critical den- 
sity; and fc is the density enhancement of the collapsing 
halo with respect to the background. Working at high red- 
shift {z = 10) has the advantage that SI ~ 1 irrespective 
of the cosmological model. Sp herical top-hat collapse then 
gives fc — 178 at virialisation (|Peacock^ll999^ . 

The temperature, Tv — 10* K, is the temperature scale 
at which hydrogen starts to coUisionally ionise. This al- 
lows for efficient atomic line cooling. Below this tempera- 
ture, hydrogen can only cool through radiative emission from 
H2 roto-vibr ational transitions, wh ich is very inefficient by 
comparison dLe Bourlot et alJll999D . Using Tv = 10* K and 
z = 10 gives M = 6.4 X lO'' Mq ~ Merit. Thus, just before 
the epoch of reionisation {z — 10), halos which reach a virial 

* In fact lMarcolini et al.l <200d) find that Merit ~ 2 X IO'^Mq as 
a result of metal cooling. This slightly lower bound doesn't affect 
our argument. 
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temperature of Tv = 10'* K and can efficiently form stars are 
those of mass Merit • 

Thirdly, t he Extended Press Schechter (EPS) scheme 
jLacev fc Colefcl993 l can be used to calculate the mean num- 
ber of halos of a given mass and redshift, M,z, which fall 
into a larger halo of given mass and redshift, A/o, zo: 
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where D = Sc [DC^)"' - ^(^o)"'] ; S = ag(Af) - crg(Afo); 
(To(Af) is the variance of the linear power spectrum at z = 
smoothed with a top-hat filter of mass M; 5c is the critical 
over-density for spherical collapse; and D{z) is the growth 
factor which depends on the cosmology^ . 

Integrating equation|5|allows us to solve for the number 
of halos, A'^, formed from the beginning of the universe up to 
some redshift z, in a given mass range, which are likely to 
have fallen into a larger halo of mass Mq = 2 x 10^^ (~ 
the mass of the MW ifwilkinson 'fc Evanj|l999ll 'l at zo = 0. 
This is what we plot in FigureQfor z = 15, 10, 5 (solid, dot- 
ted and dashed lines). The formation redshift for these ha- 
los, z, then corresponds to the epoch of reionisation, before 
which the LG dwarfs spheroidals could form stars. Over- 
plotted are the number of dSph galaxies currently observed 
around the MW (horizontal line) and the mass scale Afcrit 
(vertical line). For early reion isation z ^ 10, such as that 
favoured by the WMAP data JSoergel et alJl2003^ ■ the to- 
tal number of satellites with mass lO'^ - 10® Mq is of order 
100. This is an order of magnitude greater than the observed 
number in the Local Group. However, it is the number which 
will be accreted by a MW sized ha lo, not the number which 
survive as satellite halos at 2 = O. lMoore et al.l i2006l) have 
recently shown that only ~ 10% survive, giving the correct 
order of magnitude of satellites around the MW. A final en- 
couraging point is that the total stellar mass in the accreted 
satellites is then ~ 10®M(7), which is th e mass of the MW old 
stellar halo iBinnev fc Merrifieldll99'3) . Note that the above 
requires early reionisation. If reionisation is found to occur 
much later (2 ~ 5), this would rule out our model. Turn- 
ing this around, we tentatively suggest that the LG dwarf 
galaxies' number and distribution can constrain the epoch 
of reionisation. 

Merit ~ 10** Mq appears to be a critical mass scale at 
which gas cooling becomes efficient so that star formation 
can occur, but also at which the potential well is just deep 
enough to hold on to the remaining gas left over from star 
formation in the presence of supernovae winds. The number 
of surviving subhalos of this mass formed at 2 ~ 10 which 
fall into a MW sized halo at 2 = is consistent with the 
observed number of LG satellites; while observations of the 
LG and nearby dwarfs suggest masses of the order 10* Mq 



^ Usually it is stated that 5c = 1.69, but this is only true for 
a universe with flm = 1; where Qm is the ratio of the matter 
density of the universe to the critical density required for clo- 
sure. Here we calculate 5c and D{z) correctly (numerically) for 
a ACDM cosmology in which Qm + f^A = 1; where Q\ is the 
ratio of the dark energy density to the critical density, or the 
'cosmological constant'. The relevant equations for this are given 
in lEiscnstcin & Hu Il999|j). The cosmological parameters we use 
are given in section!^ 



(see Figure |7J . These facts, combined with the mounting 
evidence for old stellar halos in the dirr galaxies, and star 
formation histories which show pre- reionisation bursts, form 
the main motivation for investigating our model in more 
detail. 



3 HYDRO-DYNAMICAL SIMULATIONS 

We ran hydro-dynamical simulations using t he new version 
of the parallel TreeSPH code GADGET-2 dSpringel et alJ 
^01; SprinEol 2005). GADGET-2 was used in its TreePM 
mode which speeds up the calculation of long-range grav- 
itational forces considerably. The simulations were per- 
formed with periodic boundary conditions with an equal 
number of dark matter and gas particles and used the 
conservative 'entropy - formu lation' of SPH proposed by 
ISpringel fc Hernauisd (l2002ll . Radiative cooling and heat- 
ing processes were follow ed using an implementation sim- 
ilar to that of iKatz et a l. (1996) for a primordial mix of 
hydrogen and helium. We assume d a mean UV backgroun d 
produced by quasars as given bv iHaardt fc Madaul Jl99dl . 
and we switched the background on at at high redshift in 
order to have HI reionisation at 2 ~ 13. Such a treatment 
for the reionisation is quite crude: we simply stretched the 
epoch of reionisation and switched on a uniform background 
at 2 ~ 13, instead of 2 ~ 6. However, even if the details are 
incorrect (to correctly follow reionisation, one would need a 
radiative transfer code) we obtain the behaviour we expect, 
i.e. a jump in the gas temperature that brings most of the 
gas elements to T ~ 10* K. Note also that the effect of reion- 
isation only really enters at the very end of our simulations 
and serves only to inhibit further star formation. 

We used 2 x 448'^ dark matter and gas particles in a 1 
comoving Mpc//i box. The simulations were all started at 
2 — 199 and we have stored 19 redshift outputs for each 
run, mainly in the redshift range 10 < 2 < 30. The initial 
gas temperature was T = 546 K, and 40 ±2 SPH neighbours 
were used to compute physical quantities. The gravitational 
softening was set to 0.055 kpc in comoving units for all 
particles. 

Our fiducial model was a 'concordance' ACDM model 
with t7om = 0.26, f^oA = 0.74, fiob = 0.0463 and = 
72kms"^Mpc"^ n = 1 and as = O.Q*'. This model is in 
agreement with most of the observations including the pa- 
rameters inferred by the WMAP team in their first year data 
release and the recent results of the Lyman-a forest commu- 
nity (Spcracl ot al. 2003; Vicl ct al. 2004; Scliak ot al. 200^. 
The CDM transfer functions of all models have been taken 
from lEisenstein fc Hul . 



We use the common nomenclature where Honi.oA.Ob a^re the 
ratio of the density of the universe in matter, dark energy and 
baryons respectively at redshift, 2 = 0, to the critical density re- 
quired for closure; Hg is Hubble's constant at 2 = 0; n is the ini- 
tial spectral index of initial matter fluctuations (1 corresponds to 
scale-invarience) ; and erg is the amplitude normalisation of matter 
fluctuations at 2 = 0: the rms density variation of the universe 
smoothed with a top-hat filter of radius 8 Mpc {h = 0.72 
is the dimensionless Hubble parameter) at 2 = 0. Note that one 
should be cautious of this definition of trg since it implicitly as- 
sumes that linear theory can connect the initial and z = ampli- 
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Figure 2. Projected density for tlie darlc matter (left), stars (middle) and gas (rigiit), for (from top to bottom) the NWA, NWB, WWA 
and WWB simulations; second most massive halo. The gas density is for a thin slice of width O.lkpc about the x-axis. The contour bars 
show density in units of 10*' MQkpc"^ in all cases except the NW gas density, which is in units of 10'' M0kpc~^. Note that the gas is 
more extended and less concentrated than the stars and dark matter and for this reason is plotted on a different scale. All simulations 
were run with identical initial phase space distributions. This means the above halos are identical apart from their star formation and 
feedback prescriptions. The key thing to notice is that without any supernovae feedback (NWA) there is significant star formation; all 
of the most massive substructures and the central halo have bound stars. By contrast, with supernovae feedback, the NWB, WWA and 
WWB simulations are similar: there is very little star formation and significant remaining gas. The form of the feedback is not critical. 
NWB uses a heating term for the supernovae feedback; WWA uses a supernovae driven galactic wind; WWB uses both. Since NWB and 
WWA are so similar, we show only results for the NWB simulation from here on. 
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We note t hat the new three year data release by the 
WMAP team JSoergel et alj|200d) point to smaUer values 
for the spectral index and erg, while the new reionisation 
redshift is perfectly consistent with that chosen in our simu- 
lations. These new values should not have a large impact on 
the internal structure of the halos we investigate in this pa- 
per, but the lower erg will affect the statistics of such halos. 
We note that, when combined with data from the Lyman- 
Q forest, the WMAP three year releas e favours higher val- 
ues of (7s. similar to those we use here JViel et alJl2003 and 
ISeliak et al As such we felt it would be premature to 

re-run our simulations with lower as. 

We ran six different simulations in total. One with dark 
matter only was run as a control simulation and is only 
briefly mentioned. The other five investigate different star 
formation prescriptions and wind strengths. In all cases we 
started with identical initial phase space distributions. We 
used two different star formation recipes. The first (A) was 
a very simplified prescription based only on a density and 
temperature cut. All gas is turned into stars when the over- 
density, 5 > 1000 and the temperature, T < 10^ K. This 
simple prescription does not track or calculate metallicities 
for the stars and gas. The second (B) was more physically 
motivated and modelled a sub-particle multiphase medium, 
including supernovae feedback. This supernovae feedback is 
represented by a local heating term near star forming re- 
gions and does not drive a global galactic wind. This m odel 
is described in detail in ISpringel fc Hernauisd ll2003h . We 
used the 'simplified', rather than the 'explicit' mode of star 
formation. We list all the parameters for both prescriptions 
in Table [3 

In addition to each star formation prescription, we in- 
vestigated the effect of adding a supernovae driven galac- 
tic wind. The wind model we use is described in detail in 
[Soringcl & HcrriQuist (2003). We assumed that the galaxy 
mass-loss rate that goes into a wind, Mw, is proportional to 
the star formation rate itself 

Mw = r?M* , (3) 

where rj is a coefficient of order unity. Moreover, we assumed 
that the wind carries a fixed fraction x of the supernova 
energy. Equating the kinetic energy in the wind with the 
energy input by supernovae, 

^My^vl = x^snM., (4) 
we obtain the wind's initial velocity as 
/ 2X£SN 

V — V — ' 

with esN = 4 X 10*^ ergs/M©, which is the average expected 
value from the SN explosions' energy release (10®^ ergs). 
These parametrizations and the chosen values are mainly 
motivated by observations of starburst driv en galact ic winds 
in the local universe fe.g. lMartin 1999 and lott et al. 2005). 

We investigated each star formation prescription with- 
out winds (no wind A/B: NWA/B), and with winds of vary- 
ing strength (weak wind A/B: WWA/B; strong wind A; 

tude of fluctuations. On scales as large as Mpc at 2 = this 
is usually ok, but non-linearities can certai nly affect t he expected 
peculiar velocities on these scales (see e.g. |Peeill200(f) . 



SWA). The 'weak wind' (WW) and 'strong wind' (SW) runs 
had X = 0.25, 1 and rj = 2 which resulted in an average speed 
of the wind of; 221.8 and 483.6 km/s, respectively. 

Although the effects of galactic winds at 2 > 10 is not 
clear and very difficult to quantify we note that the feed- 
back prescription used here predicts a global star formation 
history and an IGM (Intergalactic Medium) metal enrich- 
ment at g = 3 that are in goo d agreement with observa- 
tions dSpringel &: Hernauisll2002'l . The effect of feedback by 
galactic winds on the IGM structures tha t surround galax- 
ies at z~2, is extremely uncertain (e.g. *Adelberger et al] 
20031. I2OO5I: Ffheu ns et al. 2002; Dcsiacques et al. 2003 
^!auchet^d]|20^f) . but we stress that the speed of the wind 
values used in the feedback simulations (although at much 
highe r redshift) are in rough agreement with local observa- 
tions JOtt et al.ll2'005ll . 

However, some caution is appropriate. There can be 
many free unconstrained parameters and, with fine tuning, 
one wonders if any set of observations could eventually be 
reproduced. To attempt to assuage these fears and give more 
confidence to the simulations, we have briefly summarised 
all of the free parameters in the model and how they have 
been constrained in Tabled the constraints are essentiall y 
the same as those outlined in lSpringel fc HernauistI i2002f) . 
The simple star formation prescription (A) is not well moti- 
vated by observations, but serves to investigate the depen- 
dence of our results on the unknown star formation recipe. 
For prescription (B), there are a number of free parameters. 
However, notice that all of these are constrained by inde- 
pendent observations except for the metallicity yield. We 
have not fine-tuned any of the parameters to achieve a par- 
ticular result. This means that our model should produce 
believable physical results, with the exception of the abso- 
lute value (not the spread) of the metallicity. The absolute 
metallicit y has been 'tuned' t o fit the metallicity of the IGM 
at 2 = 3 llTheuns et alJl2002^ . But for dwarf galaxies form- 
ing at 2 = 10, it could conceivably be quite different. We 
comment further on this in section 

It is worth briefly noting some points of contention 
about the observations we use to constrain our model pa- 
rameters. Firstly, we assume a Salpeter initial mass function 
(IMF) of stars. However, it is difficult to measure the mas- 
sive end of the stellar mass function since these stars are so 
short-lived. Yet this is the region of most interest for deter- 
mining the strength of the supernovae feedback and winds. 
Certainly the first-forming stars would appear, in theoret- 
ical models, to have an IMF biased towards the high-mass 
end i Abel et"al .1199811 . Secondly, our wind model posits that 
winds are driven during any period of star formation, even 
during relatively quiescent periods (recall that Mw = ?7M*). 
Alter natives link winds only to starbursts which are short- 
lived jMac Low fc Ferraralll99gl) . There is, then, a potential 
danger in linking observations of star bursting dwarfs with 
our more continuous wind model. However, in the simula- 
tions, the majority of the stars form over a short timescale 
(see sectional, so this should not be too great a concern. 

The galaxies were extracted from the simulation vol- 
ume using a Friends-of-Friends algorithm wit h linking length 
I = .2. We follow the iterative method of IPorciani et all 
i2002h to ensure that the extracted halos are bound. At each 
stage the total energy of each particle is calculated. Particles 
which do not appear to be bound are excluded from the po- 
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tential calculation in the next stage. The momentum centre 
of the bound group is also calculated; this means particles 
excluded at one stage may later re-enter the calculation. We 
perform our analysis using only the bound component, al- 
though our tests show the difference in profiles obtained is 
not significant. 



4 RESULTS 

In this section we present the results from our suite of hydro- 
dynamical simulations. Recall that we have to stop these 
simulations at z = 10 due to our small box size. This makes 
connecting the galaxies we form in our simulations to galax- 
ies observed at 2; = difficult. The most massive galaxies we 
form have mass ~ lO^M©. From the linear theory arguments 
given in section |21 we expect that ~ 100 of these galaxies 
will be accreted onto a galaxy the size of the Milky Way. 
Most will be accreted; ten percent will survive, as shown by 
iMoore et alJ (|200(J l. We can expect many of these survivors 
at 2 = to look very different from the galaxies we form 
stt z — 10. However, the smallest gas-poor galaxies in the 
LG at 2 = (the dSphs) are known to have very old stellar 
populations. As such, we can constrain our model by com- 
paring the stellar distributions and metallicity of our most 
massive galaxies with data from the dSphs of the LG at the 
present epoch (see sections 14.21 and 14.411 . The connection to 
the dirr galaxies will be more tentative, but we attempt this 
in section 1^31 




Figure 3. Dark matter density profiles for the most massive halo 
in the NWA, NWB and WWB simulations. Over-plotted is the 
density profile of the most massive halo in the dark matter only 
simulation (thick black line) , and a theoretical prediction for how 
this profile would respond to the adiabatic contraction of the 
baryons in the NWA simulation (thick dashed line). 



4.1 Overview 

Figure 121 shows the projected dark matter, star and gas den- 
sities for the second most massive halo in the NWA, NWB, 
WWA and WWB simulations. We do not show the results 
for the SW simulations since these are identical to the WW 
runs: the wind strength, for the plausible range of observed 
wind speeds, is not an important factor in our wind model. 
We present results for the SWA simulation alongside the 
others in section 

In the NWA simulation there is significant star forma- 
tion; all of the most massive substructures and the cen- 
tral halo have bound stars. This recovers the familiar 'over- 
cooling problem', which we discuss further in section 14.31 
There is very little gas left over and the gas density is two 
orders of magnitude lower than in all other cases. By con- 
trast, the NWB, WWA and WWB simulations are similar: 
there is very little star formation and significant remaining 
gas. The remaining gas shows bubble-like regions with over- 
densities of a factor of ~ 20. These over-densities are caused 
by heating from the star forming regions. In the case of the 
WWA model, this is a result of the supernovae driven galac- 
tic winds; for the NWB model, it is the result of heating 
from the supernovae feedback. 

The key point to take away from Figure |5| is that all of 
the simulations with feedback produce qualitatively similar 
results: suppressed and spatially extended star formation. 
In particular, the NWB and WWA simulations are almost 
identical in every respect despite their quite different star 
formation and feedback mechanisms. For this reason, we plot 
only the results for the NWB simulation from here on. We 
come to why the NWB and WWA simulations are so similar 



in section None of the simulations result in a significant 
'blow out' of gas, even in the galactic wind models. This is 
to be expected from the analytic arguments given in section 

m 

Recall that we used identical initial phase space con- 
ditions for each of the simulations. We plot identical halos 
which differ only in their star formation and feedback pre- 
scriptions. Yet the evolution of substructure in each of the 
simulations is quite different. The dark matter halo in the 
NWA simulation is much more centrally concentrated. This 
is a result of the increased merging due to the condensa- 
tion of a large mass fraction of baryons at the centre of each 
substructure halo. By contrast, the WWA and NWB simula- 
tions are much less concentrated, corresponding to less merg- 
ing. Finally, the WWB simulation (which has the strongest 
feedback - a combination of supernovae heating and galactic 
winds) shows the least merging of all. Substructure can still 
clearly be seen to the bottom left of the main halo, both in 
the dark matter and gas distributions, and below right in 
the dark matter. 

The dissipation and collaps e of gas significantly alters 
the merger history of halos (c.f. iMaccio et al]l200(j) . How- 
ever, it does not alter the dark matter density profiles, just 
the concentration. In Figure|3| we plot the dark matter den- 
sity profiles for the most massive halo in the NWA, NWB 
and WWB simulations. Over-plotted is the density profile 
of the most massive halo in the dark matter only simulation 
(thick black line), and a theoretical prediction for how this 
profile would respond to the adiabatic contraction of the 
baryons in the NWA simulation (thick dashed line) . For this 
analytic calculation, we used the usual prescription given in 
iBlumenthal et al.l Jl986l) . This assumes that the halo col- 
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S.F.P. 


Free Parameter 


Descripiion 


Constraint 


A 


5c = 1000 


Critical overdensity for star formation 


Unconstrained 


A 


Tc = 105 K 


Critical temperature for star formation 


Unconstrained 


B 


P = 0.1 


Mass fraction in stars > 8Mq 


Salpeter IMF 


B 


= 2.1Gyrs 


Star formation time 


Fit to Kennicutt ( 19891 law 


B 


P = 3Zq 


Mctallicity yield 


Based on Thcuns ct^l^ ^200^; poorly constrained 


A/B 


eSN = 10*4 J 


Supernovae energy 


Canonical value 


A/B 


7? = 2 


Wind mass loss rate 


Observed COtt et al. 2005: Martin 19991 


A/B 


''wind =0, 221.8 and 483.6 km/s 


Wind speed 


Observed (Ott et al. 2005: Martin 19991 



Table 1. Free parameters in the model and their observational constraints. Parameters are split into those which govern star formation 
(above) and those controlling the supernovae winds (below). The columns from loft to right give the star formation prescription; free 
model parameters; a brief description of the free parameter; and observational constraints on that p arameter. The metallicity is not 
constrained by independent observations (observations not of dwarf galaxies). We use a value based on lTheuns et all i2002ft . They find 
p = 3 provides the best fit to the metallicity of the IGM at ^ = 3. 




R(kpc) R(kpc) R(kpc) 




Figure 4. Projected surface density (top) and velocity dispersion (bottom) for stars in the two most massive halos (solid and dotted 
lines) in the NWA, NWB and WWB simulations. The thick and thin lines show two different projection angles which bracket the range 
of projected profiles. Over-plotted are data from two LG dSphs which br acket the observed range of brightness: Draco (crosses) and, 
the re centl y discovered, Ursa Major (squares). The data were taken from iKlevna et all ^2001^ . Iwilkinson et all <2004h . IwTUman et all 
J2003) and lKlevna et alj f2005l. Ursa Major is so faint that it has only one data point for the kinematics, but this suggests a very similar 
central velocity dispersion to Draco, despite it being an order of magnitude fainter. For Ursa Major's surface brightness distribution, we 
use the measured value for the scale length and central brightness, but assume the same distribution as in Draco - hence the lack of error 
bars for these points. Note that the data is from dSph galaxies at redshift 2 = 0, while the simulations stop at z = 10. The comparison 
is only meaningful because the dSph stellar populations are old (see start of this section for further details) . 
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lapses spherically without shell crossing, conserving mass 
and specific angular momentum and with all halo particles 
moving on circular orbits. It provides a reasonable fit in the 
central regions to the NWA simulation, but performs poorly 
over intermediate radii. This is a kno wn result. Even before 
the Blumenthal paper, lYound ijlQSffl demonstrated analyt- 
ically that the assumption that the halo particles move on 
circular orbits l eads to an over-estim ation of the contraction. 
More recently, iGnedin et al.l i200'j) have used simulations 
to demonstrate the same disparity on galaxy cluster scales. 
Here we verify that the problem persists on the scale of 
dwarf galaxies. The standard adiabatic contraction model 
should be used with caution on these scales. With super- 
novae feedback there is a lower central halo density than in 
the dark matter only run, while without feedback, the Blu- 
menthal model over-predicts the effect. This agrees well with 
previous findings in the literature (Gncdin & Zhao 2002; 
pRead fc Gilmord l2005l) . All of the dark matter halos are 
well-fit by NFW profiles llNavarro et alJll99d) : the typical 
concentration parameter and mass are c ~ 4, M ~ 10* Mq . 
Even in the presence of a strong wind, a central dark matter 
cusp persists. 

4.2 Stars viewed in projection 

Figure |1] shows the projected surface density (top) and ve- 
locity dispersion (bottom) for the stars in the two most 
massive halos (solid, dotted and dashed lines) in the NWA, 
NWB and WWB simulations. The thick and thin lines show 
two different projection angles. Over-plotted are data from 
two LG dSphs which bracket the observed range of bright- 
ness: Draco (crosses) and, the recently dis covered, Ursa 
Major (squares). The data we r e taken from [ K lgYna e t alJ 
l^2001^. I W ilkinson et all l)2004h . IWillman et alJ l^200.'^^ and 
iKlevna et al. ( 20Q5ll . Ursa Major is so faint that it has only 
one data point for the kinematics, but this suggests a very 
similar central velocity dispersion to Draco, despite it be- 
ing an order of magnitude fainter. For Ursa Major's surface 
brightness distribution, we use the measured value for the 
scale length and central brightness, but assume the same 
distribution as in Draco - hence the lack of error bars for 
these points. 

In the NWA simulation a significant amount of sub- 
structure within the main halo forms stars (see Figure I^J. 
This can clearly be seen in the surface density and projected 
velocity dispersion of the stars, both of which, for the second 
most massive halo, show a significant bump correlated with 
the position of the largest substructure halo. By contrast 
the NWB simulation, in which the substructure halos form 
few or no stars, show no such features. Such substructure is 
interesting. Several of the LG dwarfs appear to have bumps 
and wiggles in their projected velocity dispersion profiles. 
For two, the projected velocity dispersion appears to fall 
sharply to wards the edge of the light, which is dropping off 
smoothly iWilkinson et al.l2004ft . For the WWA simulation, 
we did find one halo which retained a few stars within a sub- 
halo. This can be seen in Figure |21 notice for the stars in 
the WWA simulation, there is an overdensity to the bottom 
left correlated with a dark matter subhalo. These subhalo 
stars also lead to a sharp drop in the projected velocity dis- 
persion along some lines of sight, with no associated sharp 
truncation in the stellar surface density. 



Such substructures are unlikely to be very long-lived. 
However, it is interesting to speculate that kinematic bumps 
and wiggles could correspond to late-mfalling substructure 
within the dwarfs. 

The NWA simulation halos have very high stellar sur- 
face densities and consistently falling projected velocity dis- 
persions, whereas the LG dSphs have low surface densities 
and flat projected velocity dispersions (see data points, mid- 
dle panels). This highlights the importance of supernovae 
heating and galactic winds for keeping the surface density 
of the stars low and producing galaxies which qualitatively 
resemble the dSphs of the LG. 

The NWB halo much more closely resembles the LG 
dwarfs in their stellar distribution (see middle panels). The 
agreement inside ~ 400 pc is very encouraging, especially 
given that we have made no attempt to fine tune our model 
to fit the data. However, for our simulated galaxies, the pro- 
jected surface brightness and velocity dispersion both fall 
too steeply beyond ~ 400 pc. 

The WWB model - which includes both supernovae 
heating and galactic winds - seems to solve these problems. 
The stars are more extended and the velocity dispersions 
are flat. This is exactly what we expect from gas mass loss. 
For the NWB simulation, the gas contributes significantly 
to the potential from ~ 0.4 kpc outwards. This explains why 
the stellar velocity dispersion is falling rather than flat (as 
would be expected if only the dark matter halo contributed 
to the p otential). As shown in^Read fc Gilmord ijioO^ (and 
see also iMashchenko et all l2005l) gas mass loss will cause 
the remaining stars to expand and settle into an approxi- 
mately exponential distribution - similar to that observed 
in Draco. Since the final potential is then dominated by the 
dark matter, the resulting velocity dispersions are fiat. 

However, it is important to remember that we are com- 
paring galaxies forming at z = 10 with those observed at 
2: = in the Local Group. As such, the above discrepancies 
could also be the natural result of external feedback which is 
not included in our model. The tidal field from a large nearby 
galaxy, like the MW, will heat stars beyond a characteristic 
'tidal radius' (rt), leading to more extended surface bright- 
ness profile s and flat or rising projected velocity dispersions 
beyond rt llRead_et_alJ j2006a|l . Using the analytic formulae 
for rt given'lri TR^a^^^IT l2006b ) . we find for a Draco of 
total mass 1O*M0, with a tidal radius of rt = 400 pc, that 
its orbital pericentre is ~ 20 kpc. This is co nsistent with 
measurements of the proper motion of Draco jKlevna et alJ 
I2OOII) . An alternative possibility is gas mass loss due to ram 
pressure stripping, which is likely to mimic the effect of the 
galactic wind in WWB. We return to these issues in section 
1131 

4.3 Angular momentum and star formation 

The angular momentum of a halo is typically parameterised 
by the dimensionless spin parameter given by jBuUock et al.l 
2001): 



JV 

J^ ^ \^m^,j r. X (7) 

j=0 
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Figure 5. Distribution of spin parameters (A') for the NWA, NWB and WWB simulations. The soUd, dotted and dashed histograms 
show A' for the dark matter, stars and gas. The smooth sohd line is a lognormal fit to the dark matter A' (see equation |HJ. The best-fit 
parameters are marked on each plot. 
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Figure 6. Metallicity distributions of the three most massive 
halos in the WWB simulation (thick solid, dotted and dashed 
histograms). The smooth solid lines show data for the LG dSphs 
taken from Mateo 11998). We have assumed that these distribu- 
tions are Gaussian, although detailed recent o bservations suggest 
mor e complex distrib utions are likely (see e.g. lTolstov et al 120041 
and iKoch et alJl200d) . 



where the subscript, i, denotes the particle species (dark 
matter, stars or gas), Ji is the total angular momentum in 
that species, 7?vir is the virial radius^, Mi is the mass of 
that species interior to i?vir and K-ir = \/ GMvir/ Rvir is 
the circular speed at the virial radius. This definition means 



^ The virial radius is here defined as the radius within which the 
mean density of the halo is equal to 200 times the critical density 



Figure 7. Total stellar mass as a function of dark matter mass 
within the virial radius for the ten most massive halos in the 
NWA, NWB, WWA, WWB and SWA simulations. Over-plotted 
are data for the LG dwarfs and nearby dIrrs. See text for details 
of the data compilation. 



that \[ = 1/ \/2 if all of the particles are orbiting on circular 

orbi ts at Ryji. 

iBullock et alJ (120011) demonstrated that the lognormal 
profile provides an excellent fit to the distribution of spin 
parameters in their simulations. iCoh 'n et al. (2004) showed 
that the same form provides an excellent fit to dwarf galaxy 
scale halos. The lognormal profile is given by: 



P(A') 



: exp 



ln^(A7A;,) 
2cr2 



(8) 



of the universe at . 



10. 



In Figure |K| we plot the distribution of spin parameters 
(A') for the NWA, NWB and WWB simulations. The solid. 
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dotted and dashed histograms are for the dark matter, stars 
and gas. The smooth sohd line shows the log normal fit to 
the dark matter; the best fit parameters are shown in the 
top-right of each plot. These are in excellent agreement with 
those found by IColfn et alJ ll2004) . In each case we used 
all halos with more than ~ 10 particles (in dark matter, 
stars and gas respectively), which is about 1000 halos per 
simulation. 

For the NWA simulation, we recover the well-known 
over-cooling problem. The simple star formation prescrip- 
tion we use turns all of the cold gas into stars and over- 
merging leads to a very low final stellar angular momen- 
tum. The NWB simulation alleviates this problem, but low 
angular momentum stars still remain. The WWB simula- 
tion, with the strongest feedback, prevents almost any an- 
gular momentum loss from the stars. As expected analyt- 
ically (see e.g. [Mailer & Dckcl 2002), feedback solves the 
over-cooling problem by preventing the subhalos from efh- 
cientl y forming stars. Howev er, we should be a little cau- 
tious. lKaxifmann et alJ l)200ft) have recently shown that, for 
disc galaxies, > 10® particles are required per galaxy for an- 
gular momentum loss due to spurious numerical transport 
to be at an acceptable level. Even with our very high res- 
olution, our best resolved halos have O(IO^) particles - an 
order of magnitude too low to avoid such numerical over- 
co oling. The angul ar momentum loss mechanisms identified 
bv lKaufmann et al. (20Q 6) are: viscous friction between the 
cold disc and hot halo; transfer between the spiral arms and 
bar of the disc and the dark matter halo; and angular mo- 
mentum loss from infalling cold gas clouds. In our simula- 
tions, none of these mechanisms can act. In the simulation 
without supernovae feedback (NWA), all subhalos turn their 
gas into stars. In this case, discs do not form because of an- 
gular momentum loss in mergers. With supernovae feedback 
(NWB, WWA, WWB, SWA), we prevent stars from forming 
in the smaller subhalos. Now discs do not form because the 
velocity dispersion in the gas, maintained by the feedback 
at ~ lOkm/s, is comparable to the rotational velocity. It is 
essentially a disc with a very large scale height. We can see 
this in Figure |K| where in all of the simulations with super- 
novae feedback the stars and gas retain their initial angular 
momentum, but do not significantly condense. For these sim- 
ulations, with supernovae feedback, stars form stochastically 
where the density fluctuates upwards. It is likely that such 
fluctuations are imprinted by numerical noise, but they can 
be thought of as physical. In the 'real' universe, such fluc- 
tuations are likely to exist as a result of external tides and 
internal micro-physics. 

The above explains why the feedback recipe is not crit- 
ical in determining the final morphology of the galaxy. Pro- 
vided that it prevents the subhalos from efficiently forming 
stars then it will work. Galactic winds do this by ejecting a 
proportion of the gas available for star formation; feedback 
from supernovae heating does this by heating locally cold 
gas which would oth erwise form stars . This is a known re- 
sult. For example, in lEke et alj ll200Cfl . they show that the 
over-cooling problem can be solved by simply delaying star 
formation until after most subhalos have merged into the 
main galaxy. 

Finally, note that while in the NWB and WWB simu- 
lations the stars retain significant angular momentum with 
Ay ~ 0.05, this would be very hard to detect in practise. 



None of the galaxies show measurable tip-to-tip rotation of 
greater than ~ 2km/s. This is because the stars remain ex- 
tended and do not collapse and 'spin up'. 

4.4 Metallicities 

In Figure |S] we plot the metallicity distributions of the three 
most massive halos in the WWB simulation (solid, dotted 
and dashed histograms); the NWB simulation showed very 
similar results (recall that the simple star formation pre- 
scription - A - does not produce reliable metallicity out- 
put). The smooth solid lines show data for the LG dSphs 
taken from iMated lll99^ . We have assumed that these dis- 
tributions are Gaussian, although detailed recent observa- 
tions suggest more com plex distributions are likely (see e.g. 
iTolstov et a'i]l2004l and Koch et"ai]l2006D . 

Recall from section |^ that the absolute value of the 
metallicity for our halos should be taken with caution as we 
have chosen a yield which produ ces a good fit to the metal- 
licity of the IGM at 2 = 3 (Th euns et alJl2002ll . It is not 
clear that this is the relevant value for our dwarf galaxies 
forming at z = 10. That caveat aside, our metallicities are 
systematically lower than the observed dwarfs. This is, per- 
haps, to be expected: the galaxies we form in WWB seem to 
be less massive, both in stars and dynamically, than most 
observed nearby dwarfs (see Figure 2] and Figure |7| in the 
following section). Recent observations of the faintest and 
lowest mass dwarf galaxies suggest metallicities which are 
lower than the mean of those shown in Figure 

M The faintest 

galaxy to date - Bo otes - has a mean metallicity of ~ —2.5 
l|Munoz et alJlioO^ . closer to the most massive galaxy we 
form in WWB (especially considering that the yield is so 
poorly constrained). 

The spread, which is more certain in our models, ap- 
pears to match quite well that of the nearby dwarfs. But 
in our models it is clearly non-Gaussian (recall that we as- 
sumed Gaussianity for the data). The shape of our distri- 
bution matches very well that of more deta iled modelling 
recently performed bv lMarcolini et al] (2006') . We find, sim- 
ilar to their results, an asymmetric distribution peaked at 
high me tallicity, with a tail to lower metallicities. Recent 
work bv iBellazzini et al] (|22q2) suggests such asymmetries 
may be present in the Draco and UMi dSphs, but it would 
be interesting to confirm this as a generic feature in future 
observational work. 

Our model may provide a simple solution to the 'abun- 
dance problem' for the LG dwarfs. iTolstov et al.l ll2003h 
found, from detailed spectroscopy of resolved stars, that: 1) 
the abundance^of stars in four LG dSphs is lower than that 
observed in the old stars of the MW stellar halo; and 2) the 
mean metallicity of the MW halo stars is much lower than 
that of the dwarfs: <[Fe/H]>~ —3. This presents a puz- 
zle since it suggests that the MW stellar halo cannot have 
formed from galaxies like the LG dSphs seen at the present 
epoch. In our model this is not a problem. Our progenitor 
'baryonic building blocks' (BBBs) do have the correct mean 
metallicity. Furthermore, while we do not explicitly track 
abundances in our code, we can expect their abundances 
to be high. This is because there is a link between abun- 
dance and the star formation timescale. Stars which form 
rapidly are only enriched by Type II supernovae which do 
not produce much iron; stars which form more slowly can 
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be enriched by Type la supernovae which produce more iron 
relative to other metals and there fore lower the abundance 
(see e.g. lBinnev fc Merrifield ^l998'). Our early forming BBBs 
form stars rapidly and are likely to be of high abundance. 
Over a Hubble time, as our BBBs gradually form more stars, 
they then lower their abundance and start to look more like 
the LG d Sphs. This agrees well with earlier studies in the 
literature ("Bulloc k fc Johnstonir2005l iRobertson et al.ll2005l 
and iFont et aL.200^ 



4.5 Mass to light ratios 

In Figure |7| we plot the total stellar mass as a function of 
dark matter mass within the virial radius for the ten most 
massive halos in the NWA, NWB, WWA, WWB and SWA 
simulations. We stop at the first ten halos since the smallest 
of these halos have only 10* M© in stars which corresponds 
to just ~ 10 star particles. Over-plotted are data for the LG 
dwarfs and nearby dirrs. 

It is difficult to obtain an accurate compilation for 
nearby dwarf galaxies. The dIrrs have only gas kine- 
matic measurements from which the dynamical mass can 
be derived; while the dSphs have only stellar kinemat- 
ics (since they are devoid of gas). A further compli- 
cation comes from the quality of the older versus the 
newer data. Some of the dSphs now have excellent data 
(like that over-plotted in Figure But most have older 
data which was typically taken only near the centre of 
the galaxy. It is now known that these older measure- 
ments systematically underestimate the mass of the dSph 
galaxies. Our adopted solution is to use only the latest 
data. T hese have b een compiled from Bceum & Chen ga in ii 
J2004). iBeeum fc Chengalur (2003) Begum et al. (2009)? 
iKlevna et al.l (l2005l). Willman (|2005^ ."iMunoz et alJ 



\2m^. iKlevna et alJ (l2004). I Wilkinson et al. I (l2004) 



.1 (^2001^■ IWalker et all |200fJ) . 



and IWestfall et al 



('2006")" ' Irwin fc HatzidimitrioJ " 
(2006 ). For the dSphs, we plot the masses derived from de- 
tailed mass modelling in all cases except for Ursa Major (the 
faintest dSph)^. For this dSph and for the dIrrs we assume 
an isotropic isothermal sphere model. This gives: 



M(r) 



G 



(9) 



where G is the gravitational constant, the measured mean 
projected velocity dispersion squared, r = 1 kpc ~ _Rvir is 
the radius enclosing the mass and / is a small correction 
factor. / = 1.4 is chosen such that we would recover the 
correct mass for the Draco dSph (obtained from detailed 
distribution function modelling) using equation]^ 

There are two key things to take away from Figure |7| 
Firstly, notice how rapidly the stellar mass falls as the halo 



® The abundance is the mass of a given 'metal' (element heavier 
than Helium) relative to iron. 

^ Ursa Major has only one velocity dispersion measurement. It 
could be reasonably argued that this makes it as unreliable as the 
older data. However, unlike the older data, this one data point 
lies ~ 250 pc away from the centre of the galaxy. This makes it 
more likely to be representative of the mean velocity dispersion 
than older measurements. 



mass is reduced in all of the simulations with internal feed- 
back from supernovae. The details of the star formation and 
feedback recipes are not important for determining the fi- 
nal stellar masses. In all cases, the stellar mass falls two 
orders of magnitude for a decrease in halo mass of little over 
30%. This was expected from the analytic arguments pre- 
sented in section|5| Secondly, notice that the LG dwarfs and 
nearby dIrrs seem to show a cut-off in their dynamical mass 
at ~ 5 X IO'^Mq, but a wide range in stellar masses. The 
data are well bracketed by our NWA and other simulations. 

The dSph galaxies are consistent with our most mas- 
sive halos with all of the gas removed; while the dirr lie 
close to the NWA model suggesting that, over a Hubble 
time, they have managed to turn most of their initial gas 
into stars. That this process is a slow one is important. 
Firstly, because star formation is inefficient, it allows gas 
to be slowly removed. Ram pressure strippin g is one such 
mech anism which could achieve this (see e.g. iMaver et all 
|200^. It has the advantage that galaxies which lie close to 
the MW or M31 will be more rapidly stripped than those 
which lie further away, naturally reproducing the distance 
morphology relation (see section . Secondly, recall that 
in the NWA simulation, the stellar density was much higher 
than that observed in the LG dwarfs and too much substruc- 
ture formed stars. This occurred because the star formation 
in the NWA simulation was too efficient. Even for the most 
isolated dIrr which eventually turns all of its gas into stars, 
supernova feedback and reionisation have a role to play in 
keeping the star formation efficiency low and spatially ex- 
tended. This prevents substructure forming stars and keeps 
the stellar surface density and rotational velocity low, con- 
sistent with observations. 

One caveat is worth mentioning. In all cases, it is our 
most massive simulated galaxies which resemble the least 
massive LG dSphs. For the WWB simulation, none of our 
simulated galaxies achieve enough stellar mass to be consis- 
tent with observations (this can also be seen in Figure UJ. 
This suggests that either our feedback prescription is too 
strong, or the observed dwarfs formed from the merger of 
several of our BBBs, with some continued star formation 
over a Hubble time. It is not possible to address which of 
these is correct without running larger-box simulations of 
comparable resolution and continuing down to lower red- 
shift. We leave this for future work. 



5 DISCUSSION 

5.1 A connection to galaxies in clusters? 

iTrentham fc HodgkinI (|20o3) have recently discovered a new 
class of very low surface brightness (VLSB) galaxies in the 
Virgo cluster. These galaxies are extremely faint, spheroidal 
in morphology and extended over several kpc. So far in this 
paper, we have talked only about the low density environ- 
ment of the LG. In clusters, the picture could be quite dif- 
ferent. While little is known about these VLSB galaxies and 
follow-up observations are currently underway, an intrigu- 
ing possibility is that these too are naked stellar halos. Like 
the dSph galaxies of the LG, they may have lost their gas 
faster than they could turn it into stars, leaving them with 
just the extended, low rotational velocity, metal poor, old 
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stellar halo component. If true, further observations should 
confirm that these VLSB galaxies have all of the properties 
of old stellar halos. They are likely more massive than their 
LG cousins, however. 



5.2 The formation of globular clusters 

We have presented the case for a 'baryonic building block' 
for galaxy formation. These gas-rich building blocks have 
total mass Merit ~ 1O*M0 and form ~ lO^M© in stars be- 
fore reionisation. However, we can see in the universe that 
stars form on smaller scales than this in star clusters and 
globular clusters. In part this is a natural result of fragmen- 
tation: in our model, these star clusters and globular clusters 
form within such baryonic building blocks. However, the pic- 
ture cannot be quite this simple for two reasons. Firstly, we 
require that the phase space density of stars in our build- 
ing blocks is low, yet in globular clusters and massive star 
clusters it is very high. Secondly, we have outlined the im- 
portance of supernovae winds in suppressing star formation 
on mass scales ^ IO^Mq. How then can star clusters with 
stell ar mass ^ lO^ Mp, ever form? 

'Kr oupa fc Boilv (,2002.1 provide some plausible answers 
to these problems. On the scale of star clusters, the star for- 
mation timescales are shorter than the lifetime of a typical 
0-star. This explains why gas expulsion from supernovae is 
not likely to be of great importance for star cluster forma- 
tion (though it becomes very important for driving galactic 
scale winds and heating surrounding gas, as we have shown 
here). Instead, stellar winds and ionising flux from O-stars 
are the most important forms of internal feedback for star 
clusters forming with masses < 1O*M0. Above lO^Mo, su- 
pernovae become important for star clusters and can allow 
them to self enrich. Such massive clusters are still protected 
from destruction by supernovae winds, however, by their 
high density (Morgan fc Lak c 1989). 

It is not clear how the very high gas densities re- 
quired to form globular clusters can be achieved. Observa- 
tional evidence suggests a link to the star formation rate. In 
interacting systems (like galaxy mergers) and strong star 
bursts, the globular cluster and massive star cluster for- 
mation rate is much la rge r than in more quiesce nt systems 
JWhitmore et al.ll993l and lZepf fc Ashmanll99^ . This view 
is als o supported by recent theoretical work (see e.g. lLi et alJ 
I2OO4I) . It is clear that accurately modelling such effects is be- 
yond the scope of this work. 



6 CONCLUSIONS 

Using analytic arguments and a suite of very high resolu- 
tion (~ 10'^ M0 per particle) cosmological hydro-dynamical 
simulations, we have argued that high redshift, z ~ 10, 
M ~ 10^ M0 halos, form the smallest 'baryonic building 
block' (EBB) for galaxy formation. These halos are just mas- 
sive enough to efficiently form stars through atomic line cool- 
ing and to hold onto their gas in the presence of supernovae 
winds and reionisation. These combined effects, in particu- 
lar that of the supernovae feedback, create a sharp transi- 
tion: over the mass range 3 — 10 x lO^M©, the BBBs drop 
two orders of magnitude in stellar mass. Below ~ 2 X 1O^M0 



galaxies will be dark with almost no stars and no gas. Above 
this scale is the smallest unit of galaxy formation: the BBB. 

We have shown that supernovae feedback is important 
for these smallest galaxies. Not because it ejects the gas, 
but because it keeps the gas hot and extended. We find that 
the details of such feedback is not critical. Whether imple- 
mented as a heating term for gas surrounding star forming 
regions, or as a galaxy-wide wind, the results are similar. 
However, the combination of supernovae heating and super- 
novae driven galactic winds gives the best agreement with 
observations. Such feedback works by reducing the star for- 
mation efficiency in subhaloes. This keeps the surface den- 
sity and rotational velocity of the stars which do form low. 
The smallest observed galaxies in the Local Group have very 
low surface brightness, as does the Milky Way old stellar 
halo. Without such supernovae feedback our model cannot 
reproduce these properties. Efficient cooling of the super- 
novae ejecta is prevented by the ionising background from 
reionisation. 

We connected these BBBs to galaxies observed at 2 = 
using a mixture of linear theory arguments and results 
from other studies in the literature. In a galaxy the size 
of the Milky Way, O(IOO) such building blocks will be ac- 
creted. iMoore et alj i200(J) have recently shown that, of 
these, ~ 10% survive to form the lowest mass Local Group 
dwarf galaxies. The remainder form the bulk of the Milky 
Way old stellar halo. The survivors will slowly form stars and 
lose gas over a Hubble time. Since neither reionisation nor 
supernovae winds actually eject gas from these galaxies, we 
require some other mechanism for this. Ram pressure strip- 
ping is a likely candidate. In this case, those BBBs on benign 
orbits which keep them far away from the Milky Way or An- 
dromeda manage to retain their gas and slowly form stars - 
these become the lowest mass dwarf irregular galaxies (dirr) ; 
those on more severe orbits lose their gas faster than they 
can form stars and become the dwarf spheroidals (dSphs). 
Both galaxy types become more metal rich and of lower 
abundance than the progenitor building blocks due to their 
more extended star formation. In this picture, the dSphs are 
'naked stellar halos', while the dirrs have both an old, metal 
poor, stellar halo and a younger, more centrally concentrated 
population of stars. There is incr easing observational evi- 
dence that this is indeed the case jMinniti fc ZiilstrallT99^ . 
iMinniti et al]l2003l and lAparicio et alJbOOoT 

We have shown that the stars in our high redshift BBBs 
resemble in those in the faintest dSphs in the Local Group at 
the present epoch. This does not prove, how ever, that all of 
the Lo cal Group dSphs formed in this wav. IKravtsov et alJ 
I2OO4) and iMaver et alJ (|200^ have suggested that some 
formed from quite different mechanisms and were much more 
massive in the past than at present. 

There should be many galaxies in the Local Group 
with surface brightness an order or two orders of magni- 
tude fainter than those already found. We predict that these 
galaxies will have very similar total mass (2 — 10 x 10^ M©) 
to those satellites already discovered. 

Finally, we have commented almost exclusively on the 
low-density Local Group, rather than the high-density clus- 
ter environment. This is mainly because the BBBs we form 
in our simulations are so small and faint that, at the mo- 
ment, we can only hope to observe their z = Q counterparts 
in the very nearby universe. We expect, however, many such 
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BBBs to survive to the present epoch, even in cluster en- 
vironments. The continued accretion of such BBBs over a 
Hubble time should lead to an old rnctal-poor stellar halo 
being a ubiquitous feature of all large galaxies. However, it is 
important not to think of these BBBs as the building block 
for galaxy formation. We stress that they are the smallest 
building block and the major contributor to old stellar halos, 
not the bricks from which all galaxies are made. 
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